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(54) Title: EXAMINATION OF BREAST TISSUE USING TIME- RESOLVED SPECTROSCOPY 
(57) Abstract 

A method and a system (4) for breast tissue examination includes a time- 
resolved spectroscopy apparatus (20; 20A; 20B), a support (8; 9; II; 12; 13) 
with an input port (14) and an output port (16) separated by a selected distance 
is positioned relative to the examined breast. Locations of the input and output 
ports are selected to examine a tissue region of the breast. A light source (32; 34; 
60) generates pulses of electromagnetic radiation of a selected wavelength in the 
visible or infrared range. The pulses of duration on the order of a nanosecond or 
less are introduced into the breast tissue at the input port (14) and detected over 
rime at the detection port (16). Signals corresponding to photons of detected 
modified pulses are accumulated over time. Values of the scattering coefficient 
or the absorption coefficient of the *T«m;™H breast tissue are calculated based 
on the shape of the modified pulses. The examined breast tissue is characterized 
based on the values of the scattering coefficient or the absorption coefficient. 
Absorbing or floorescing contrast agents may be introduced into the examined 
tissue. The system may be adapted for use with x-ray mammography, needle 
localization procedure or MRI mammography. 
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EXAMINATION OF BREAST TISSUE 
USING TIME-RESOLVED SPE CTROSCOPY 

Background of the Invention 
5 The invention features a time-resolved 

spectroscopic method and apparatus for breast tissue 
examination. 

Breast cancer is among the most common and the 
most feared malignancies in women. It has an 

10 unpredictable course, the treatment is frequently 
physically and emotionally draining and the risk of 
metastatic spread persists for many years. Due to its 
high occurrence rate, routine breast cancer screening, 
which includes physical examination and x-ray 

15 mammography, plays an important role in current health 
care. X-ray mammography can detect over 90% of all 
masses and increases the 10-year survival rate to about 
95% for patients with cancers solely detected by 
• mammography. Although the modern mammography uses a low- 

20 dose of x-rays, it still involves some small risk of 

inducing cancers by the radiation. Other tests, such as 
magnetic resonance imaging (MRI) and gadolinium enhanced 
MRI, have been used successfully for detection of breast 
tumors and may be used routinely for screening in the 

25 future. 

After a small suspicious mass is detected in the 
breast non-invasively, excisional biopsy is usually 
performed to exclude or diagnose malignancy- The biopsy 
specimen is removed under local anesthesia and is used 
30 for histopathological diagnosis. The statistics show 
that in about 75% of the excisional biopsies, the 
biopsied tissue is diagnosed to be benign. Thus, a 
majority of patients undergoes this unpleasant and costly 
procedure unnecessarily. Furthermore, it has been 
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suggested that the excisional biopsy may cause spreading 
of the malignant tumor cells. 

Therefore, a non-invasive, relatively inexpensive 
technique that can detect and characterize breast tumors 
5 may find its place in today's health care, alone or in 
conjunction with the above-mentioned techniques. 

Summary of the Invention 
The invention features a system and a method for 
breast tissue examination using time-resolved 
10 spectroscopy. 

In general, in one aspect, the method includes 
the following steps. A support that includes an input 
port and an output port separated by a selected distance 
is positioned relative to the examined breast. Locations 
15 of the input and output ports are selected to examine a 
tissue region of the breast. Light pulses of a selected 
wavelength and duration less than a nanosecond are 
introduced into the breast tissue at the input port and 
detected over time at the detection port. Signals 
20 corresponding to photons of detected modified pulses are 
accumulated over the arrival time of detected photons. 
Values of a scattering coefficient or an absorption 
coefficient of the examined breast tissue are calculated 
based on the shape of the modified pulses. The examined 
25 breast tissue is characterized based on the values of the 
scattering coefficient or the absorption coefficient. 

In general, in another aspect, the method 
includes the following steps. A support that includes an 
input port and an output port separated by a selected 
30 distance is positioned relative to the examined breast. 
Locations of the input and output ports are selected to 
examine a tissue region of the breast. Light pulses of a 
selected wavelength and duration less than a nanosecond 
are introduced into the breast tissue at the input port 
35 and detected over time at the detection port. Signals 
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corresponding to photons of detected modified pulses are 
integrated over at least two selected time intervals 
separately spaced over the arrival time of the modified 
pulses. A value of an absorption coefficient of the 
5 examined breast tissue is calculated based on the shape 
of the modified pulses. The examined breast tissue is 
characterized based on the value of the absorption 
coefficient. 

In this aspect, the method may include further 

10 steps. The detected photons are integrated over other 
selected time intervals separately spaced over the 
arrival time of the modified pulses. Time dependence of 
the light intensity is determined based on the number of 
photons integrated over each time interval, and a value 

15 of a scattering coefficient of the examined breast tissue 
is determined. The examined breast tissue is 
characterized based on the value of the scattering 
coefficient. 

Preferred methods use the above-described steps 
20 and additional steps as follows. 

The input port and the output port are moved to a 
different location to examine another tissue region of 
the breast. Values of the scattering coefficient or 
absorption coefficient are again determined by repeating 
25 the above-described steps for the newly selected tissue 
region. The tissue region is characterized using the 
additional values of the scattering coefficient or the 
absorption coefficient. 

The above-described steps are performed over 
30 several tissue regions to examine the entire breast. 

The characterizing step includes comparing the 
calculated values of the scattering or absorption 
coefficient with selected values of scattering or 
absorption coefficient, respectively. 
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The selected values of the scattering and 
absorption coefficient correspond to normal breast 
tissue, normal contralateral breast tissue or series of 
homogenous breast tumors. 
5 The characterizing step includes comparing the 

calculated values of the scattering coefficient or the 
absorption coefficient with selected values of the 
scattering coefficient or the absorption coefficient, 
respectively . 

10 When the characterizing step reveals that the 

examined tissue includes abnormal tissue further the 
steps are performed. Another location of the input port 
and the output port is selected to define a new tissue 
region proximate to the region having abnormal tissue. 

15 The values of the scattering coefficient or the 

absorption coefficient of the newly selected tissue 
region are determined by applying the corresponding 
steps. Abnormal breast tissue is localized by comparing 
values of the scattering coefficient or the absorption 

20 coefficient of different selected tissue regions. The 
type of the abnormal tissue may be determined by 
comparing values of the scattering coefficient or the 
absorption coefficient of the localized tissue to values 
of the scattering coefficient or the absorption 

25 coefficient corresponding to selected tissue masses. 

The tissue masses include one of the following: 
carcinoma, fibroadenoma or fibrocystic tissue. 

The size and location of the abnormal tissue 
region is determined. 

30 When the characterizing step reveals that the 

examined tissue includes abnormal tissue further the 
steps are performed. A contrast agent exhibiting known 
optical properties at the selected wavelength is injected 
into the blood stream of the subject. Another location 

35 of the input port and the output port is selected to 
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define a new tissue region proximate to the region having 
abnormal tissue. The values of the scattering 
coefficient or the absorption coefficient of the newly 
selected tissue region are determined. The abnormal 
5 breast tissue is localized by comparing values of the 
scattering coefficient or the absorption coefficient of 
different selected tissue regions. 

The type of the abnormal tissue may be determined by 
comparing values of the scattering coefficient or the 

10 absorption coefficient of the localized tissue to values 
of the scattering coefficient or the absorption 
coefficient corresponding to selected tissue masses 
comprising the contrast agent. 

When the characterizing step reveals that the 

15 examined tissue includes abnormal tissue further the 
steps are performed. A contrast agent exhibiting known 
optical properties at the selected wavelength is injected 
into the abnormal tissue. A location of the input port 
and the output port is selected. The values of the 

20 scattering coefficient or the absorption coefficient of 
the newly selected tissue region are determined. The 
abnormal breast tissue is localized by comparing values 
of the scattering coefficient or the absorption 
coefficient of different selected tissue regions. 

25 The type of the abnormal tissue may be determined 

by comparing values of the scattering coefficient or the 
absorption coefficient of the localized tissue to values 
of the scattering coefficient or the absorption 
coefficient corresponding to selected tissue masses 

30 comprising the contrast agent. 

The contrast agent is a fluorescing material or 
absorbing material. The contrast agent is preferentially 
absorbed by the tissue mass. 
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The above described steps are performed in 
conjunction with x-ray mammography, MRI mammography or a 
needle localization procedure. 

In another aspect, the invention features a system 
5 for performing the above-described method. 

Brief Description of the Drawing 
Fig. 1 depicts diagrammatical ly a time-resolved 
spectroscopic system for breast tissue examination. 
Figs. 1A, IB, 1C and ID depict different 
10 embodiments an optical fiber support for breast tissue 
examination. 

Fig. 2 depicts diagrammatically a single photon 
counting TRS apparatus arranged for breast tissue 
examination. 

15 Fig. 3 depicts diagrammatically a TRS boxcar 

apparatus arranged for breast tissue examination. 

Fig. 3A shows a timing diagram of the apparatus of 

Fig. 3. 

Fig. 3B shows a typical time resolved spectrum 
20 collected by the apparatus- of Fig. 3. 

Fig. 4 depicts diagrammatically examination of 
breast tissue using a fluorescing contrast agent. 

Figs. 4A and 4B depict diagrammatically 
examination of breast tissue using MRI and time-resolved 
25 spectroscopy. 

Figs. 5A, 5B, 5C, 5D, 5E and 5F display values of 
the absorption coefficient and the scattering coefficient 
of normal breast tissue measured at different locations 
of the right breast and of the left breast. 
30 Figs. 6A and €B display values of the absorption 

coefficient and the scattering coefficient, respectively, 
of normal breast tissue for women of different ethnic 
background. 
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Description of the Preferred E mbodiments 
Fig. 1 depicts a breast tissue examination system 
4 placed on a human breast 5 for breast tissue 
examination. The system includes an optical fiber 
5 support 8 with multiple input ports 14 and multiple 

output ports 16. Support 8 is placed around breast 5 so 
that input ports 14 and output ports 16 define 
irradiation locations and detection locations on the skin 
of breast 5, respectively. Connected to selected input 
10 and output ports are optical fibers 15 and 17, 

respectively. System 4 uses either a single photon 
tissue resolved apparatus 2 OA or a time resolved 
apparatus 2 OB using boxcar type integration. 

Referring also to Figs. 1A, IB, 1C and ID, system 
15 4 uses different types of the optical fiber supports 
designed to introduce and detect photons at selected 
locations and thus shape the optical field. The optical 
fiber supports are made of flexible or rigid materials 
and are shaped to accommodate breasts of different 
20 volumes. Furthermore, the inside surface of the 

supports may include material of known scattering and 
absorptive properties. The material is selected to 
either return back to the breast tissue photons escaping 
through the skin (i.e., a low absorber and high 
25 scatterer) or provide additional paths for the escaping 
photons to the detector (i.e., the material has 
substantially the same optical properties as normal 
breast tissue) . The supports are designed for use with a 
time-resolved spectrophotometer (TRS) alone or in 
30 conjunction with x-ray mammography, MRI or a needle 

localization procedure. Specifically, fiberoptic support 
9 shown in Fig. 1A includes three sets of the input and 
detection ports labeled 10a, 10b, and 10c. Sets 10a and 
10c are used to measure control data and set 10b is used 
35 to examine a suspected mass 7. Furthermore, support 9 
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enables precise characterization of the distances between 
the three sets, between input ports 14 and detection 
ports 16 (p) and from the chest wall 6 to each set 
Supports 11 and 12, shown in Figs. IB, 1C and ID, are 
5 used with x-ray mammography and needle localization 
procedure, respectively, and their functions are 
described below. 

Referring to Fig. 2, a dual wavelength, time 
correlated single photon counting TRS apparatus 2 OA is 
10 connected to support 13 positioned on breast 5. Pulsed 
laser diodes 32 and 34 (model PLP-10 made by Hamamatsu, 
Japan) , are driven by a 5 mW pulser 36 connected to a 100 
MHz pulse generator 37, and generate light pulses on the 
order of 500 psec or less. The light from laser diodes 
15 32 and 34 is electro-mechanically time shared using a 60 
Hz vibrating mirror 33 and is coupled to one end of 
optical fiber 15. Optical fiber 15, which has about 200 
Mm diameter, alternatively conducts pulses of 754 nm and 
810 nm light to input port 14. The introduced photons 
20 migrate in the examined breast tissue and some of them 
arrive at output port 16. Optical fiber 17 collects 
photons of the modified pulses from an area of about 10 
mm 2 and transmits them to a PMT detector 40. 

The output of PMT 40 is connected to a wide band 
25 amplifier 42 with appropriate roll-off to give good pulse 
shape and optimal signal to noise ratio, output signals 
from amplifier 42 are sent to a high/ low level 
discriminator 44, which is a pulse amplitude 
discriminator with the threshold for pulse acceptance set 
30 to a constant fraction of the peak amplitude of the 
pulse. Next, the discriminator pulses are sent to a 
time-to-amplitude converter (TAG) 46. TAG 46 produces an 
output pulse with an amplitude proportional to the time 
difference between the start and stop pulses received 
35 from pulser 36. The TAC pulses (47) are routed by a 
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switch 48 to either a multichannel analyzer (MCA) 50 or 
an MCA 52. Switch 48 operates at 60 Hz and is 
synchronized with mirror 33. The photon emission, 
detection cycle is repeated at a frequency on the order 
5 of 10 MHz. 

Each MCA collects only a single photon for each light 
pulse introduced to the tissue. Each MCA acquires and 
sums photons of only one designated wavelength and stores 
the corresponding pulse of a shape that depends on 
10 properties of the examined tissue. The pulses are 

preferably accumulated over about 2 to 3 minutes so that 
at least 10 s counts are collected at the maximum of the 
pulse shape. The detected pulse shape is analyzed by a 
computer 56. Computer 56 is connected to pulse generator 
15 37 and MCAs 50 and 52 via an interface module 54 and is 
adapted to control the entire operation of the system. 

Alternatively, TRS apparatus 20 represents a 
boxcar TRS apparatus 20B, as shown in Fig. 3. A pulsed 
laser diode 60 is driven by a 5 mW pulser 62 connected to 
20 a 100 MHz pulse generator 64. Laser diode 60 generates a 
train of 100 ps light pulses of 754 nm wavelength coupled 
to optical input fiber 15. The light pulses are 
introduced to breast tissue at input port 14. The 
introduced photons migrate in the examined tissue and a 
25 portion of them arrives at an output port 16. In the 
migration process, the input pulse has been modified by 
the scattering and absorptive properties of the examined 
tissue. Photons arriving at detection port 16 are 
transmitted by optical fiber 17 to a detector 66, (for 
30 example, Hamamatsu photomultipliers R928, R1517, MCP 
R1712, R1892) . 

The output of detector 66 is amplified in a wide 
band preamplifier/ impedance changer 67 and coupled to a 
boxcar integrator 68. Integrator 68 activated by a pulse 
35 gate 73 collects all arriving photons over a 
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predetermined time interval 75, as shown in Fig. 3A. The 
integrator output (78) is sent to computer interface 
module 80 and computer 82. Computer 82 stores the total 
number of counts detected during the collection interval 
5 of integrator 68. 

Integrator 68 includes a trigger 70, which is 
triggered by a signal 63 from pulser 62. Trigger 70 
activates a delay gate 72 which, in turn, starts counting 
of all detected photons during the time interval 

10 specified by a gate width 74. Output from a gate width 
norma lizer 76 is an analog signal or a digital signal 
representing all photons that arrived at detection port 
16 during the preselected gate width interval (75) . A 
suitable integrator is a boxcar SR 250 manufactured by 

15 Stanford Research Systems. 

Depending on the application, computer 82 sets the 
delay time (71) of delay gate 72 and the gate width time 
(75) of gate width circuit 74. Gate width normalizer 76 
adjusts the width of the integration time depending on 

20 the detected signal level. The gate width may be 

increased logarithmically for signals at t»t max , wherein 
the detected number of photons decreases exponentially; 
this increases the signal-to-noise ratio. Furthermore, 
computer 82 can scan the integration gate widths over the 

25 whole time profile of the detected pulse. By scanning 
the delay times (71) and appropriately adjusting the gate 
widths (75) , the computer collects data corresponding to 
the entire detected pulse. Subsequently, computer 82 
calculates the shape (85) of the detected pulse and 

30 stores the time dependent light intensity profile I(t). 
The pulse shape, I (t), detected either by 
apparatus 2 OA or apparatus 2 0B possesses information 
about the scattering and absorptive properties of the 
examined breast tissue and is used to determine the 

35 scattering and absorption coefficients. Referring to 
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Fig. 3B, a test measurement was performed on apparatus 
20A and found that due to a somewhat slow response time, 
the detector broadens the reflectance profile, as seen on 
spectrum 86. Thus the experimental spectra (87) are 
5 deconvoluted to separate the instrumental response from 
the profile dispersion due to the diffusion. The 
deconvolution yields about 6% increase in the value of p a 
and about 23% decrease in the value of M 8 - 

The examined tissue region is defined by the 

10 distribution of photon pathlengths forming an optical 
field in the tissue. The size and shape of the optical 
field is a function of the input-output port separation 
(p) as well as the optical properties of the tissue 
(i.e., absorption coefficient, scattering 

15 coefficient, p e , and the mean cosine of anisotropic 

scattering, g) . The general diffusion equation is used 
to describe the photon migration in tissue, as analyzed 
by E.M. Sevick, B. Chance, J. Leigh, S. Nioka, and M. 
Maris in Analytical Biochemistry 125, 330 (1991), which 

20 is incorporated by reference as if fully set forth 
herein. The diffusion equation is solved for the 
intensity of detected light in the reflectance geometry, 
R(p,t), or the transmittance geometry T(p,d,t) . In the 
reflectance geometry, in a semi-infinite media with the 

25 separation of the input and output ports on the order of 
centimeters, the absorption coefficient is a function of 
the reflectance spectrum as follows: 

— logJ?(p,t) = — - u.c + P 2 (l) 

For t -♦ » the absorption coefficient n a is determined as 
follows : 
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lim-^log^tp, t) = - u fl c (2) 

wherein p is the separation between input and detection 
ports and c is speed of light in the medium. However, it 
is difficult to measure the at t»t max because in this 
region the data show substantial noise. Thus, to measure 
5 M a at t>>t max' requires determination of the pulse shape 
at a high number of counts. 

If the approximation of infinite time is not 
valid, Eg. 1 can be rewritten to obtain M a as follows: 

^=-^log^(p. t ) (3, 



The value for D can either be the average value obtained 
10 from numerical simulations or a value specific to the 
type of tissue being measured. 

The effective scattering coefficient (l-g) 'ti Q is 
determined as follows: 

wherein t,^ is the delay time at which the detected 
15 reflectance time profile (R(p,t) e I(t>) reaches maximum. 
After detecting the pulse shape corresponding to the 
examined tissue the computer calculates the absorption 
coefficient (/i a ) , and the scattering coefficient (M 8 ) • 
The absorption coefficient is quantified by evaluating 
20 the decaying slope of the detected pulse, using Eqs. 2 or 
3. The effective scattering coefficient, (l-g) 'P e , is 
determined from Eq. 4. 

The breast screening procedure starts by selecting 
a support with appropriate arrangement of input ports 14 
25 and output ports 16. The absorptive and scattering 
properties of the tissue are measured for one set of 
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ports and then the optical field is transferred by using 
another set of ports. The entire breast is examined by 
selecting sequentially different ports. In the 
reflection geometry, the optical field can be represented 
5 by a three dimensional "banana-shaped" distribution 
pattern or, in the transmission geometry, a "cigar- 
shaped" distribution pattern. In the "banana-shaped" 
pattern, the shallow boundary is due to the escape of 
photons that reach the air-scatter er interface while the 

10 deeper boundary is due to attenuation of long path 

photons by the absorbers. The penetration depth of the 
photons is about one half of the port separation (p) . 
During the screening procedure, the computer calculates 
/i a and p B for the entire breast and compares, the measured 

15 values with threshold values of /i a and /i B of normal tissue 
or series of /* a and ji 8 values of different homogeneous 
tumor types. As is shown in Figs. 5A through 5F, from 
one person to another there is some variation in *i a and fi a 
for normal tissue, but only a very small variation 

20 between the left breast and the right breast of the same 
person. Cancerous tissue, which is usually highly 
perfused, exhibits higher values of jt a and p B than fibrous 
tissue. Normal tissue, which has a relatively high 
amount of fat, exhibits the lowest values of M a ^s* 

25 Alternatively, instead of calculating p a and p B , 

the system can calculate an average pathlength of the 
migrating photons. From the detected and deconvoluted 
photon intensity profile, R(t) , a mean pathlength of the 
distribution of pathlengths <L> is determined as follows: 




(5) 



WO 95/02987 



PCT/US94/07984 



- 14 - 

wherein c is the speed of light in vacuum and n » 1.36 is 
the average refractive index of tissue. 

If a breast tumor is outside of the optical field, 
it does not alter the banana-shaped distribution of 
5 path lengths. As the optical field is moved closer to the 
breast tumor, which is a strongly absorbing mass, the 
photons that have migrated the farthest distance from the 
input and detection ports are eliminated by the 
absorption process. Since photons with the longest 

10 pathlengths are absorbed, the system detects reduction in 
the average path length. When the optical field is moved 
even closer to the mass, some of the detected photons now 
migrate around the mass without being absorbed; this is 
detected as lengthening of the distribution of 

15 pathlengths. Thus, the average pathlength measurement 
can reveal location of the breast mass. 

In the screening process, the breast tissue is 
characterized by several tissue variables which may be 
used alone or in combination. System 20 measures the 

20 absorption coefficient, the scattering coefficient, blood 
volume or tissue oxygenation using one or more selected 
wavelengths. The wavelengths are sensitive to naturally 
occurring pigments or contrast agents that may be 
preferentially absorbed by the diseased tissue. Suitable 

25 naturally occurring pigments are, for example, hemoglobin 
(Hb) or oxyhemoglobin (Hb0 2 ) sensitive to 754 run and 816 
nm, respectively. 

Alternatively, suitable color dyes, such as 
cardio-green or indocyinin-green, may be injected to the 

30 blood system alone or bound to a vehicle such as a 
gadolinium contrast agent, which is preferentially 
absorbed by tumors in the first five to ten minutes. An 
appropriate wavelength is selected for the color dye, for 
example, cardio-green exhibits maximum absorption at 805 

35 nm. 
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The computer can create "maps" of the breast by- 
mapping the spacial variation of the measured values for 
p at fi Qt blood volume or hemoglobin saturation. The 
resolution is enhanced when several tissue variables are 
5 mapped. The blood volume is measured using a pair of 
contrabestic wavelengths (e.g., 754 nm and 816 nm) or the 
isobestic wavelength (i.e., 805 nm) . The hemoglobin 
saturation (Y) is measured at two wavelengths (e.g., 754 
nm and 816 nm) and is calculated by taking the ratio of 
10 absorption coefficients at these wavelengths and then 
using the following equation: 

ruioo%)= ^ (6) 

wherein the coefficients are determined from the 
extinction values of hemoglobin at 754 nm and 816 nm that 
are e m> = 0.38 cm" 1 mM" 1 , = 0.18 ctBThtiJC 1 , 

15 respectively, and the difference in extinction 

coefficients between oxyhemoglobin and hemoglobin that 
are Ae^q.j^ = 0.025 cm" 1 mM" 1 and Ae^Q.^ = 0.03 cm" 1 mM" 1 , 
respectively. 

In another preferred embodiment, TRS apparatus 20 

20 is used in combination with x-ray mammography. The 
combined procedure is performed if a suspected mass is 
detected by the above-described optical method, x-ray 
mammography or other screening method. 

Referring to Pig. IB, breast 5 is compressed in 

25 either a horizontal or vertical position between an x-ray 
film case 90 with the x-ray film and a support 11 with 
input ports and output ports located on a grid. An x-ray 
mammogram is taken to determine location of suspected 
mass 7 relative to the grid. Suitable input port 14 and 

30 output port 16 are selected so that the introduced 
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optical field 92 encompasses mass 7. Then, TRS apparatus 
2 OA or 20B is used to measure /t a , /t a , blood volume or 
oxygen concentration of the examined tissue using the 
above-described techniques. The measured values are 
5 again compared to the values corresponding to normal 
tissue or different types of diseased tissue to 
characterize the mass. If an unequivocal result is 
obtained, exploratory excisional biopsy is not needed. 

In another preferred embodiment, TRS apparatus 20 

10 is used in combination with the needle localization 

procedure. The needle localization procedure locates the 
mass that is then examined by system 4. Furthermore, a 
needle used in the needle localization procedure may 
introduce an optical fiber directly to mass 7. 

15 Referring to Figs. 1C and ID, breast 5 is 

compressed between x-ray film case 90 and a support 12 
with input ports and output ports located on a grid and a 
centrally located opening for needle 94 or 98. One or 
more x-ray mammograms are taken to determine the location 

20 of suspected mass 7 relative to the grid and to the 

needle opening. The needle is inserted into the breast 
and the needle tip is positioned in the center of mass 7. 
Additional x-ray mammograms may be taken to verify or 
adjust the position of the needle. 

25 As shown in Fig. 1C, if the needle is used only to 

localize mass 7, input port 14 and output port 16 are 
selected so that their separation is equal or larger than 
two times the depth of mass 7. This separation assures 
that the introduced optical field 96 encompasses mass 7. 

30 After needle 94 is positioned, a tiny wire is inserted 
into the mass and left there for marking purposes. TRS 
apparatus 20 measures fi & , ft B , blood volume or oxygen 
concentration of the examined tissue using the above- 
described techniques. The measured values are again 

35 compared to the values corresponding to normal tissue and 
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different types of diseased tissue, and the tissue of 
mass 7 is characterized. 

As shown in Fig. ID, needle 98 locates an end 99 
of optical input fiber 15 directly inside mass 7. The 
5 optical fiber with a diameter of about 100 /im or less is 
threaded inside needle 98. End 99 is slightly extended 
from the needle so that the introduced photons are 
directly coupled to mass 7. The location of detection 
port 16 defines optical field 100. In this arrangement, 

10 all detected photons migrate in the targeted tissue; this 
increases the relative amount of the targeted tissue 
being examined and thus increases the resolution of the 
system. To compare tissue of mass 7 with normal tissue, 
the same geometry of the input port and the detection 

15 port is used to measure the optical properties of the 

contralateral breast. Alternatively, in the same breast, 
needle 98 is moved outside of mass 7 so that the 
positions of end 99 and detection port 16 define an 
optical field completely removed from mass 7. To 

20 characterize the mass, the values of /* a and m b measured 
for mass 7 are compared either to the values of normal 
tissue measured in the same arrangement or to values of 
different typ es of diseased tissue. 

In another embodiment, detection contrast is 

25 enhanced using fluorescing contrast agents. A tumor is 
permeated with the fluorescing contrast agent that has a 
decay time of less than 1 nsec, and the labeled tumor is 
then examined using a TRS device 20. Suitable agents 
emit fluorescing radiation in the range of 800 nm to 1000 

30 nm, for example, carbocyanene dyes or porphorin 
compounds. 

Referring to Fig. 4, the fluorescing contrast 
agent is injected to the blood system alone or bound to a 
vehicle, such as a gadolinium contrast agent, which is 
35 initially preferentially absorbed by a breast tumor. 
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Alternatively, the fluorescing agent is injected directly 
into the tumor. TRS device 20 generates 150 psec light 
pulses introduced into breast 5 at input port 14. The 
introduced photons of selected excitation wavelength 
5 reach tumor 7 and excite a fluorescing radiation 110, 
which propagates in all directions from tumor 7. The 
photons of fluorescing radiation 110 migrate in the 
examined breast . tissue and some of them arrive at output 
port 16. Output port 16 has an interference filter that 

10 passes only photons at, the wavelength of fluorescing 

radiation 110. Optical fiber 17 collects the transmitted 
photons, which are delivered to the PMT detector. System 
4 may detect fluorescing radiation 110 at several output 
ports at the same time or move the ports to different 

15 positions on the fiber optic support. 

TRS device 20 detects the pulses of fluorescing 
radiation 110; the shape of these pulses depends on the 
decay time of the fluorescing agent and optical 
properties of both the tumor tissue and the normal breast 

20 tissue. 

In another embodiment, referring to Figs. 4A and 
4B, TRS apparatus 20 is used in combination with HRI. 
MRI examined the breast with or without using a rare 
earth contrast agent. A network of surface coils 112 is 

25 cooperatively arranged with a fiberoptic support 114, 
which is constructed for use with MRI. The network of 
coils 112 and support 114 are appropriately located 
around the examined breast. At the same time as the MRI 
data are collected, TRS apparatus 20 collects the optical 

30 data. If an abnormal mass is detected, the MRI 
identifies the size and location of the mass. The 
optical data are then used to characterize the mass. 
Optical contrast agents may be used alone or in 
combination with the rare earth contrast agents, as 

35 described above. 
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EXPERIMENTS: 

Preliminary experiments were conducted under a 
pre-approved protocol and after receiving informed 
consent of women with normal breast tissue and of women 
5 having a mass detected in their breast. 

The examination of normal breast tissue was 
performed at several different locations of the right and 
left breast of the same woman. Referring to Figs. 2 and 
3, letters RR, LR, LR and LL denote the right and left 

10 breast and the right and left breast side where the input 
and detection ports were located, respectively. Figs. 5A 
and 5B summarize the absorption coefficient (M a ) anci tne 
scattering coefficient (fi e '= (1-g) *M B ) , respectively, 
measured on the right and left breast at a separation p = 

15 6 cm. The values of pt a and p B ' for the right breast are 
identical to the values for the left breast within the 
measurement error. Figs. 5C and 5D summarize fi B ' , 
respectively, for the tissue on the left side and the 
right side of each breast, and Figs. 5E and 5F summarize 

20 /i a , /i B ', respectively, for the tissue located at different 
distances from the chest wall. The data shown in Figs. 
5A through 5F confirm that there is no significant 
difference in the optical properties measured over the 
entire breast. 

25 The examination of normal breast tissue also 

measured differences in the optical properties 
corresponding to the volume of the breast, type of the 
breast tissue, the age of the woman and ethnic 
background. Based on the decreasing X-ray background 

30 absorption, the breast tissue was categorized as "dense", 
"fatty", or a "mixture" of the two tissue types. The 
values of fi & and p B * of "fatty" tissue are lower than for 
"dense" tissue, which has a higher fibrous content. 
Since the shape of the breast varies, it was difficult to 
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categorize precisely the breast volume. For the volume 
measurement, the breast was stabilized on a plate and the 
length was measured from the chest to the end of the 
breast. The width and the thickness were measured at 
5 approximately l cm from the chest. Tissue of a large 
volume breast exhibits lower values of fi a and /x, 8 ' than 
that of a small volume breast. The same trend is 
observed for women above the age 50 when compared to 
women below 50. Figs. 6A and 6B display values of fi A and 

10 /i a ' , respectively, for Caucasians and African-Americans. 
Normal breast tissue of Caucasians and African-Americans 
exhibits substantially the same optical properties except 
for the values of fi a ' measured at the 4 cm separation. 
Since the skin forms a higher relative percentage of the 

15 examined tissue at a smaller separation of the input and 
output ports, the lower values of /a b ' may be due to a 
lower scattering coefficient of the skin with more 
pigment . 

In all measurements, a smaller separation of the 
20 input and output ports yielded larger values of fi a and fi e ' 
than a larger separation of the ports. This differences 
can be explained by violation of the semi-infinite 
boundary conditions at the smaller separations, i.e., a 
larger escape of photons through the tissue surface 
25 before they are collected by optical fiber 17. 

Furthermore, this dependence exists since the slope of 
the photon decay was measured not sufficiently far from 
the peak of the reflectance data as expressed in Eqs. 1, 
2 and 3. This problem arises due to a low photon count 
30 of approximately 10,000 counts at the peak. Thus the 
measured data have a low signal to noise ratio and a 
reliable reflectance data can not be taken at t»t max . 
The corresponding error in the absorption coefficient, 
EfPrtHabg, is determined using Eq. 7. 
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B(p , t)/ ^ = At- (7, 

The error of scattering coefficient, E(p,t, /J a ) | BCt , 
arises due to the error in the absorption coefficient. 
The corresponding error is determined using Eq. 8. 

Eip, t,\i a )/ set «- i- t4E(p,t)/ ai)S C 2 t l Lc]-S(p, t)/^ (8) 
3p' 



The preliminary values of M a and /i B ' corrected for the 
5 error using Egs. 7. and 8 are shown in Table 1. 



TABLE 1 


Separation 
(cm) 


Mean 


Ma 

Error 
Adjusted 


Mean 
Ma' 


Mean 
tmax (ns) 


Error 
Adjusted 


4 


0.029 


0.020 


16.2 


1.55 


14.1 


) 5 


0.023 


0.021 


11.4 


1.9 


10.9 


6 


0.019 


0.022 


9.5 


2.3 


10.1 


7 


0.017 


0.021. 


7.9 


2.7 


8.7 



The corrected values of mean n a for different 
separations, p, are substantially the same, but the 

15 corrected values of mean /i e still are p dependent 
although their spread is reduced considerably. 

The examination of a breast with abnormal tissue 
was performed substantially the same way as the above- 
described examination of normal breast tissue. The 

20 breast tissue was first characterized by x-ray 

mammography and the size and location of a mass was 
determined. The examined breast was compressed between 
x-ray film case 90 and a support 12, as shown in Fig. 1C. 
Input port 14 and output port 16 were selected so that 

25 mass 7 was located in optical field 96. . 

The values of pt a and ft 8 ' measured around tumor 7 
(using input output port set 10b of Fig. 1A) were 
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compared to control data measured on the same breast 
(using input output port sets 10a and 10c of Fig. 1A) . 
The measured data were also correlated with pathology 
information on abnormalities in the examined breasts. 
5 The abnormalities were divided into the following three 
categories: fibrocystic, Fibroadenoma, and Carcinoma. 
Furthermore, these three categories are subdivided 
according to the size of the tumors as follows: smaller 
in diameter than 1 cm, and equal or larger in the 
10 diameter than 1 cm. Preliminary data measured on over 
fifty patients show an increase in the values of both fi a 
and m 8 ' when compared with normal tissue but statistical 
significance has not been demonstrated. 

Other embodiments are within the following claims: 
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1. A system for in vivo examination of breast 
tissue of a female subject, comprising: 

a light source constructed to generate pulses of 
electromagnetic radiation of a selected wavelength in the 
5 visible or infra-red range, said pulses having duration 
on the order of a nanosecond or less, 

a support, positionable relative to said breast, 
comprising 

an input port and an output port separated by 
10 a selected distance, said distance specifying a volume of 
the examined tissue, said input port, constructed to 
define an irradiation location of said breast, adapted to 
introduce into the breast tissue said generated pulses of 
light of said selected wavelength, said output port 
15 constructed to define a detection location of said 
breast, 

optical material, at least partially 

surrounding said tissue, for limiting escape of photons 

outside said tissue, 
20 a detector constructed to detect over time photons 

of modified pulses that have migrated in said breast 

tissue to said output port and produce corresponding 

electrical signals, 

a photon counter, connected to said detector, 
25 constructed to accumulate over time said electrical 

signals and determine a shape of said modified pulses, 

and 

a processor adapted to calculate values of 
scattering coefficients or absorption coefficients of the 
30 examined breast tissue based on said shape of said 

modified pulses and determine a physiological property of 
the breast tissue based on said coefficients. 

2. A system for in vivo examination of breast 
tissue of a female subject, comprising: 
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a light source constructed to generate pulses of 
light of a selected wavelength in the visible or infra- 
red range, said pulses having duration on the order of a 
nanosecond or less, 
5 a support, positionable relative to said breast, 

comprising 

an input port and an output port separated by a 
selected distance, said input port, constructed to define 
an irradiation location of said breast, adapted to 
10 introduce into the breast tissue said generated pulses of 
light of said selected wavelength, said output port 
constructed to define a detection location of said 
breast, 

optical material, at least partially 

15 surrounding said tissue, for limiting escape of photons 
outside said tissue, 

a detector constructed to detect over time, 
at said output port, photons of said modified pulses that 
have migrated in said breast tissue to said output port 

20 and produce corresponding electrical signals, 

gated integrator and an integrator timing control 
adapted to integrate said photons over at least two 
selected time intervals separately spaced over the 
arrival time of said modified pulses; and 

25 a processor adapted to calculate values of 

absorption coefficients of the examined breast tissue 
based on the number of photons integrated over each time 
interval and determine a physiological property of the 
examined breast tissue based on said coefficients* 

30 3. The system of claim 2 wherein said gated 

integrator, said integrator timing control, and said 
processor are further adapted to determine a delay time 
^ t max) between the time a pulse is introduced and the time 
the intensity of detected into the tissue corresponding 

35 modified pulse has a maximum value, said processor being 
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further adapted to determine the effective scattering 
coefficient (1-g) *M S of ^ & examined tissue. 

4. A system for in vivo examination of breast 
tissue of a female subject, comprising: 
5 a light source constructed to generate pulses of 

electromagnetic radiation of a selected wavelength in the 
visible or infra-red range, said pulses having duration 
on the order of a nanosecond or less, 

a support, positionable relative to said breast, 
10 comprising 

an input port and an output port separated by 
a selected distance, said distance specifying a volume of 
the examined tissue, said input port constructed to 
define an irradiation location of said breast, said 
15 output port constructed to define a detection location of 
said breast, 

optical material, at least partially 

surrounding said tissue, for limiting escape of photons 

outside said tissue, 
20 a detector constructed to detect over time photons 

of modified pulses that have migrated in said breast 

tissue to said output port and produce corresponding 

electrical signals, 

a photon counter, connected to said detector, 
25 constructed to accumulate over time said electrical 

signals corresponding to a shape of said modified pulses, 

and 

a processor, connected to said photon counter, 
adapted to calculate values of a scattering coefficient 
30 or an absorption coefficient of the examined breast 
tissue by comparing said detected radiation to said 
introduced radiation, and determine a physiological 
property of the breast tissue based on said coefficients. 

5. The system of claim 1, 2, 3 or 4 wherein said 
35 selected wavelength is sensitive to a contrast agent of 
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known optical properties, said contrast agent being 
injectable into tissue of said subject, and said detected 
photons of said pulses have been modified by said 
contrast agent and said tissue while migrating in said 
5 tissue to said output port. 

6. A system for in vivo examination of tissue of 
a subject, comprising: 

a light source constructed to generate pulses of 
electromagnetic radiation of a selected wavelength in the 
10 visible or infra-red range, said pulses having duration 
on the order of a nanosecond or less, 

said selected wavelength being sensitive to a 
contrast agent of known optical properties, said contrast 
agent being injectable into tissue of said subject, 
15 an input port, constructed to define an 

irradiation location of the examined tissue, adapted to 
introduce into the tissue said generated pulses of light 
of said selected wavelength, 

an output port constructed to define a detection 
20 location of the examined tissue, 

a detector constructed to detect over time photons 
of pulses modified by said contrast agent and said tissue 
while migrating in said tissue to said output port, said 
detector producing electrical signals corresponding to 
25 said detected photons, 

a photon counter, connected to said detector, 
constructed to accumulate over time said electrical 
signals corresponding to a shape of said modified pulses, 
and 

30 a processor, connected to said photon counter, 

adapted to calculate values of a scattering coefficient 
or an absorption coefficient of the examined tissue by 
comparing said detected radiation to said introduced 
radiation, and determine a physiological property of the 

35 tissue based on said coefficients. 
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7. A system for in vivo examination of tissue of 
a subject, comprising: 

a light source constructed to generate pulses of 
electromagnetic radiation of a selected wavelength in the 
5 visible or infra-red range, said pulses having duration 
on the order of a nanosecond or less, 

said selected wavelength being sensitive to a 
contrast agent of known optical properties, said contrast 
agent being injectable into tissue of said subject, 
10 an input port, constructed to define an 

irradiation location of the examined tissue, adapted to 
introduce into the tissue said generated pulses of light 
of said selected wavelength, 

an output port constructed to define a detection 
15 location of the examined tissue, 

a detector constructed to detect over time photons 
of pulses modified by said contrast agent and said tissue 
while migrating in said tissue to said output port, said 
detector producing electrical signals corresponding to 
20 said detected photons, 

gated integrator and an integrator timing control 
are adapted to integrate said photons over at least two 
selected time intervals separately spaced over the 
arrival time of said modified pulses; and 
25 a processor adapted to calculate values of 

absorption coefficients of the examined breast tissue 
based on the number of photons integrated over each time 
interval and determine a physiological property of the 
examined breast tissue based on said coefficients. 
30 8. The system of claim 1, 2, 3, 4, 5, 6 or 7 

wherein said support is constructed to enable accurate 
selection of said irradiation locations and said 
detection locations relative to said breast, each said 
selection forming a defined relative geometry of said 
35 irradiation and detection locations. 
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9. The system of claim 8 wherein said support is 
constructed for cooperative use with x-ray mammography. 

10. The system of claim 8 wherein said support is 
constructed for cooperative use with needle localization 

5 procedure. 

11. The system of claim 10 wherein said support 
further comprises an opening adapted to accomodate a 
needle used in said needle localization procedure. 

12. The system of claim 42 wherein said support 
10 is constructed for cooperative use with MRI mammography. 

13. The system of claim 1, 2, 3, 4, 5, 6, 7 or 8 
wherein said processor is further adapted to compare said 
calculated values of the scattering coefficient or the 
absorption coefficient to selected values of the 

15 scattering coefficient or the absorption coefficient, 
respectively. 

14. The system of claim 5, 6, 7 or 8 wherein said 
contrast agent is preferentially absorbed in an abnormal 
tissue. 

20 15. The system of claim 5, 6, 7 or 8 wherein said 

contrast agent is fluorescing when irradiated and said 
detector is adapted to detect fluorescent radiation. 

16. The system of claim 15 wherein said detector 
includes a filter adapted to pass only said fluorescent 
radiation. 

17. The system of claim 5, 6, 7 or 8 wherein said 
contrast agent is also magnetically active and suitable 
for MRI examination of breast tissue. 

18. The system of claim 17 wherein said contrast 
agent includes a rare earth contrast agent. 

19. The system of claim 18 wherein said contrast 
agent is preferentially absorbed by abnormal tissue. 

20. The system of claim 1, 2, 3, 4, 5, 6, 7 or 8 
wherein said contrast agent is preferentially absorbed in 
an abnormal tissue region. 
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21. A method for examination of breast tissue of a 
female subject using pulses of light of a selected 
wavelength, said method comprising the steps of: 

(a) providing a support, positionable relative to 
the examined breast, comprising 

an input port and an output port separated 
by a selected distance, said input port constructed to 
define an irradiation location of said breast, said 
output port constructed to define a detection location of 
said breast, 

optical material, at least partially 
surrounding said tissue, for limiting escape of photons 
outside said tissue, 

(b) selecting locations of said input and output 
ports to examine a tissue region of the breast, 

(c) introducing into the breast tissue, at said 
input port, pulses of light of a selected wavelength in 
the visible or infra-red range, said pulses having 
duration on the order of a nanosecond or less, 

(d) detecting over time, at said detection port, 
photons of modified pulses that have migrated in the 
breast tissue, 

(e) accumulating, over arrival time of said 
detected photons, electrical signals corresponding to 
said detected photons of said modified pulses, 

(f) calculating values of a scattering coefficient 
or an absorption coefficient of the examined breast 
tissue by comparing said detected radiation to said 
introduced radiation, and 

(g) characterizing the examined breast tissue 
based on said values of the scattering coefficient or the 
absorption coefficient. 

22, A method for examination of breast tissue of 
a female subject using pulses of light of a selected 
wavelength, said method comprising the steps of: 
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(a) providing a support, positionable relative to 
the examined breast, comprising 

an input port and an output port separated by 
a selected distance, said input port constructed to 
define an irradiation location of said breast, said 
output port constructed to define a detection location of 
said breast, 

optical material, at least partially 
surrounding said tissue, for limiting escape of photons 
outside said tissue, 

(b) selecting locations of said input and output 
ports to examine a tissue region of the breast, 

(c) introducing into the breast tissue, at said 
input port, pulses of light of a selected wavelength in 
the visible or infra-red range, said pulses having 
duration on the order of a nanosecond or less, 

(d) detecting over time, at said detection port, 
photons of modified pulses that have migrated in the 
breast tissue, 

(e) integrating said photons over at least two 
selected time Intervals separately spaced over the 
arrival time of said modified pulses, 

(f ) calculating a value of an absorption 
coefficient of the examined breast tissue based on the 
number of photons integrated over each time interval, and 

(g) characterizing the examined breast tissue 
based on said value of the absorption coefficient. 

23. The method of claim 22 wherein said (e) and 
(f) steps further comprise: 

integrating said photons over other selected time 
intervals separately spaced over the arrival time of said 
modified pulses, 

determining the time delay (t max ) between the time 
a pulse is introduced into the tissue and the time the 
intensity of corresponding modified pulse has a maximum 
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value, based on the number of photons integrated over 
said time intervals, and 

calculating a value of the effective scattering 
coefficient of the examined breast tissue. 

24. A method for examination of tissue of a 
subject using pulses of light* of a selected wavelength, 
said method comprising the steps of: 

providing an input port and an output port 
separated by a selected distance, said input port 
constructed to define an irradiation location of said 
tissue, said output port constructed to define a 
detection location of said tissue, 

selecting locations of said input and output ports 
to examine a tissue region of said subject, 

injecting into the tissue a contrast agent of 
known optical properties at said selected wavelength, 

introducing into the tissue, at said input port, 
pulses of light of a selected wavelength in the visible 
or infra-red range, said pulses having duration on the 
order of a nanosecond or less, 

detecting over time, at said detection port, 
photons of pulses modified by said contrast agent and 
said tissue while migrating in said tissue to said output 
port, 

accumulating, over arrival time of said detected 
photons, electrical signals corresponding to said 
detected photons of said modified pulses, 

calculating values of a scattering coefficient or 
an absorption coefficient of the examined tissue by 
comparing said detected radiation to said introduced 
radiation, and 

characterizing the examined tissue based on said 
values of the scattering coefficient or the absorption 
coefficient. 
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25. A method for examination of tissue of a 
subject using pulses of light of a selected wavelength, 
said method comprising the steps of: 

(a) providing an input port and an output port 
separated by a selected distance, said input port 
constructed to define an irradiation location of said 
tissue, said output port constructed to define a 
detection location of said tissue, 

(b) selecting locations of said input and output 
ports to examine a tissue region of said subject, 

(c) injecting into the tissue a contrast agent of 
known optical properties at said selected wavelength, 

(d) introducing into the tissue, at said input 
port, pulses of light of a selected wavelength in the 
visible or infra-red range, said pulses having duration 
on the order of a nanosecond or less, 

(e) detecting over time, at said detection port, 
photons of modified pulses that have migrated in the 
tissue, 

(f ) integrating said photons over at least two 
selected time intervals separately spaced over the 
arrival time of said modified pulses, 

(g) calculating a value of an absorption 
coefficient of the examined tissue based on the number of 
photons integrated over each time interval, and 

(h) characterizing the examined tissue based on 
said value of the absorption coefficient. 

26. The method of claim 25 wherein said (e) and 
(f) steps further comprise: 

integrating said photons over other selected time 
intervals separately spaced over the arrival time of said 
modified pulses, 

determining the time delay (t,^) between the time 
a pulse is introduced into the tissue and the time the 
intensity of the corresponding modified pulse has a 
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maximum value, based on the number of photons integrated 
over said time intervals, and 

calculating a value of the effective scattering 
coefficient of the examined tissue. 

27. The method of claim 21, 22, 23, 24, 25 or 26 
further comprising the steps of: 

moving said input port and said output port to a 
different location to examine another tissue region of 
the breast, 

repeating steps (c) through (f) to determine 
values of the scattering or absorption coefficient of 
said newly selected tissue region, and 

characterizing the examined breast tissue based on 
said values of the scattering coefficient or the 
absorption coefficient. 

28. The method of claim 21, 22 or 23 further 
comprising, the step of: 

injecting into the blood stream of said subject a 
contrast agent of known optical properties at said 
selected wavelength, wherein said detected photons of 
said pulses have been modified by said contrast agent and 
said tissue while migrating in said tissue to said output 
port. 

29. The method of claim 24, 25, 26, 27 or 28 
further comprising the step of comparing said calculated 
values of the scattering coefficient or the absorption 
coefficient to selected values of the scattering 
coefficient or the absorption coefficient, respectively. 

30. The method of claim 29 wherein said contrast 
agent is preferentially absorbed by abnormal tissue. 

31. The method of claim 29 wherein said contrast 
agent is fluorescing when irradiated, and said detecting 
step includes detecting fluorescent radiation. 

32. The method of claim 31 wherein said detecting 
step includes detecting only said fluorescent radiation. 
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33. The method of claim 29 wherein said contrast 
agent is also magnetically active and said method is used 
with MRI. 

34. The method of claim 33 wherein said contrast 
agent includes a rare earth contrast agent. 

35. The method of claim 29 wherein said selected 
values of the scattering and absorption coefficient 
correspond to normal breast tissue. 

36. The method of claim 29 wherein said selected 
values of the scattering and absorption coefficient 
correspond to normal contralateral breast tissue. 

37. The method of claim 29 wherein said selected 
values of the scattering and absorption coefficient 
correspond to series of homogenous breast tumors. 
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